FTUAM 04/18 
IFT-UAM/CSIC-04-43 
hep-ph/0407352 
July 2004 



^ ■ 

g . A comparison between direct and indirect 

(N ; dark matter search 

Oh' 
(D 
00 



m 
(N 



Y. Mambrini and C. Munoz 



^ ■ Departamento de Fisica Teorica C-XI and Institute de Fisica Teorica C-XVI, 

CN '. Universidad Autonoma de Madrid, Cantoblanco, 28049 Madrid, Spain. 

in 

m 

o 
o 



Abstract 



• We carry out a comparison between different kinds of experiments for the detection 

of neutralino dark matter. In particular, we compare the theoretical predictions for 
Qh! the direct detection of these particles, through their elastic scattering on target nuclei 

through nuclear recoils, with those for their indirect detection through the gamma 
rays produced by their annihilation in the galactic center. For the comparison we pay 
special attention to points which will be testable by very sensitive experiments, such 
X ; as e.g. EDELWEISS II and GENIUS for direct detection, and GLAST for indirect 

■ detection. We perform the analysis in the framework of a general SUGRA theory, 

allowing the presence of non-universal soft scalar and gaugino masses in the parameter 
space. As a conclusion, we find that indirect detection experiments will be able to test 
larger regions of the parameter space than direct detection ones. 
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1 Introduction 



A weakly interacting massive particle (WIMP) is one of the most interesting candi- 
dates for dark matter in the Universe. The main reason being that WIMPs can be 
left over from the Big Bang in the right amount to account for the observed matter 
density. Let us recall that the observations of galactic rotation curves [1], cluster of 
galaxies, and large scale flows [2], imply 0.1 < iloMh^ ^ 0.3. In addition, the recent 
data obtained by the WMAP satellite [3] confirm that dark matter must be present 
with 0.094 < QdmIi'^ < 0.129. 

Many direct detection experiments have been carried out around the world in order 
to detect WIMPs [4] . These particles cluster gravitationally in our galaxy, and therefore 
there must be a flux of them on the Earth. This flux can be detected by observing the 
elastic scattering of WIMPs on target nuclei through nuclear recoils. In fact, one of 
the experiments, the DAMA collaboration, reported data favouring the existence of a 
WIMP signal [5]. Taking into account uncertainties on the halo model, it was claimed 
that the preferred range of the WIMP-nucleon cross section IS (T 10-6 - 10-5 pb 
for a WIMP mass smaller than 500 — 900 GeV [6, 5]. Unlike this spectacular result 
other collaborations such as CDMS [7] and EDELWEISS [8], claim to have excluded 
important regions of the DAMA parameter space. 

This controversy may be solve in the future by these and other projected experi- 
ments. For example, GEDEON [9] will be able to explore positively a WIMP-nucleon 
cross section cr > 3 x 10"^ pb. Similar cross sections will also be tested by EDELWEIIS 
II. CDMS Soudan (an expansion of the CDMS experiment in the Soudan mine), will be 
able to test a > 2 x 10-^ pb. In fact, while this paper was in preparation, the Soudan 
experiment reported its first results [10], and these are even stronger than those of 
EDELWEISS. Let us finally mention that the most sensitive detector will be GENIUS 
[11], which will be able to test a WIMP-nucleon cross section as low as o" ~ 10-^ pb. 

On the other hand, there are also promising methods for the indirect detection of 
WIMPs through the analysis of their annihilation products [12]. One of them consists 
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of detecting the gamma rays produced by these annihilations in the galactic halo. For 
this, one uses atmospheric Chcrenkov telescopes or space-based gamma-ray detectors. 
Similarly to the case of direct detection, one of the space-based experiments, the En- 
ergetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma-Ray 
Observatory, detected a signal [13] that apparently cannot be explained with the usual 
gamma-ray background. In particular, after five years of mapping the gamma-ray sky 
up to an energy of about 20 GeV, there are evidences of a source of gamma rays above 
about 1 GeV in the galactic center region with a value for the integrated flux of 
about 10~^ cm~^ s~^. It is worth noticing, however, that alternative explanations for 
this result have been proposed modifying conveniently the standard theory of galactic 
gamma ray [14]. 

As in the case of direct detection, projected experiments might clarify the situation. 
For example, a larger sensitivity will be provided by the Gamma-ray Large Area Space 
Telescope (GLAST) [15], which is scheduled for launch in 2006. In addition, energies 
as high as 300 GeV can be measured testing models with heavy dark matter more 
easily than EGRET. Even if alternative explanations arise for the events observed by 
EGRET, GLAST will be able to detect a smaller flux of gamma rays from dark matter 

~ lO"^'^ cm~^ s~^. Besides, ground-based atmospheric Cherenkov telescopes are 
complementary to space-based detectors. For example, the High Energy Stereoscopic 
System (HESS) experiment [16], which has begun operations with four telescopes"*^, has 
a large effective area of 10^ m^ (versus 1 m^ for GLAST) allowing the detection of very 
high energy fluxes > 1 TeV. On the contrary, the energy threshold is larger, > 60 
GeV. 

Given the experimental situation described above, concerning the detection of dark 
matter, it is natural to wonder how big the scattering (annihilation) cross section for 
its direct (indirect) detection can be. Obviously, this analysis is crucial in order to 
know the possibility of detecting dark matter in the experiments. However, the answer 
to this question depends on the particular WIMP considered. The leading candidate 
in this class of particles is the lightest neutralino, Xi, which appears in supersymmetric 
(SUSY) extensions of the standard model [18]. 

For the case of direct detection, the analysis of neutralino dark matter has been 
carried out with the following results. In the usual minimal supergravity (mSUGRA) 
scenario, where the soft terms of the minimal supersymmetric standard model (MSSM) 
are assumed to be universal at the unification scale, Mqut ~ 2 x 10^^ GeV, and ra- 
diative electroweak symmetry breaking is imposed, the cross section turns out to be 
small, <^^0-p ^ 3 X 10~® pb [4]. Clearly, in this case, current direct-detection experi- 
ments are not sufficient and more sensitive detectors producing further data, such as 
e.g. GENIUS, are needed. However, this result can be modified by taking into ac- 
count possible departures from the mSUGRA scenario. In particular, a more general 
situation in the context of SUGRA than universality, the presence of non-universal 
soft scalar and gaugino masses, may increase the cross section significantly (for recent 

^ After this paper was sent for publication, the first data from HESS appeared [17]. Although they observe 
gamma rays in a particular direction, these would correspond to a neutralino heavier than 12 TcV, a range 
disfavoured by particle physics where supersymmetric particules should not be heavier than about 1 TeV. 
Thus, an astrophysical source is a more likely explanation for these data. 
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works see e.g. [19, 20, 21] and references therein). 

For the indirect detection of dark matter through gamma rays, similar theoretical 
analyses have been carried out for the case of mSUGRA [12, 22, 23, 24, 25, 26, 27]. 
Neutralino annihilations in the galactic center region seem not to be able to explain 
the excess of gamma-ray fluxes detected by EGRET. In fact, only small regions of the 
parameter space of mSUGRA may be compatible with the sensitivity of GLAST and 
HESS, for standard dark matter density profiles such as the Navarro, Frenk and White 
(NFW) one. However, as in the direct detection case, when non-universal scalar and 
gaugino masses are allowed, the annihilation cross section may increase significantly, 
producing larger gamma-ray fluxes [27]^. 

Given the above theoretical results, it would be interesting to carry out a comparison 
of the sensitivity of both kinds of experiments, namely direct versus indirect detection 
of dark matter. This is precisely the aim of this paper. In Section 2 we will review very 
briefly the above results concerning the direct and indirect detection of neutralino dark 
matter in a general SUGRA theory. Taking into account this discussion, in Section 3 
we will compare the sensitivity of both kinds of experiments in the parameter space of 
SUGRA. In particular, we will pay special attention to points which will be testable by 
very sensitive experiments, such as e.g. EDELWEISS II, GENIUS, GLAST. Finally, 
the conclusions are left for Section 4. 

2 SUGRA predictions for the detection of dark 
matter 

In this section we will briefly review the SUGRA scenario and its impact on the direct 
and indirect detection of neutralino dark matter. Let us first recall that in mSUGRA, 

where the soft terms are assumed to be universal, and radiative electrowcak symmetry 
breaking is imposed, one has only four free parameters defined at the GUT scale: the 
soft scalar mass m, the soft gaugino mass M, the soft trilineax coupling A, and tan /?. 
In addition, the sign of the Higgsino mass parameter, fj,, remains also undetermined by 
the minimization of the Higgs potential. 

Using these parameters the neutralino-proton cross section has been analysed ex- 
haustively in the literature. The result is that the scalar cross section is small, and 
only future detectors such as GENIUS will be able to cover some regions of the pa- 
rameter space (for a recent analysis, see e.g. [19]). Similar results are obtained when 
gamma-ray fluxes from the galactic center are studied in the framework of mSUGRA. 
In particular, the annihilation cross section of neutralinos is small, and we need future 
detectors such as GLAST to analysed at least some regions. 

However, as discussed in detail in [19, 20, 21], in the context of direct detection, 
and in [27] for indirect detection, both cross sections mentioned above can be increased 
significantly when the structure of mSUGRA soft terms is abandoned"^ Therefore 

^For analyses in the context of the efiMSSM scenario, where the parameters are defined directly at the 
electroweak scale, see e.g. the recent works [28, 29, 30]. 

■^In this sense, it is worth noticing that such a non-universal structure can be recovered in the low-energy 
limit of some phenomenologically appealing string scenarios. This is the case, for example, of D-brane 
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larger regions of the parameter space can be covered by the experiments. Following 
those works, there are several interesting possibilities when non-universal scalars and 
gauginos are considered in a general SUGRA theory. 

On the one hand, parameterising the non-universalities in the Higgs sector, at the 
GUT scale, as follows: 

mji^ = m\l + 6i) , mjj^ = m\l + 62) , (2.1) 

one has three representative cases: 



a) Si = 


, S2 


= 1 , 




b)Si = -1 


, S2 


= 0, 




c)Si = -1 


, S2 


= 1 . 


(2.2) 



On the other hand, parameterising the non-universalities in the gaugino sector as fol- 
lows: 

Ml = M , M2 = M(l + 62) , M3 = M(l + 5'^) , (2.3) 

where Mi 2,3 are the bino, wino and gluino masses, respectively, one has the represen- 
tative case: 

d) S'^ = -0.4 , S'2 = 0. (2.4) 

Of course, as discussed in detail in [21], the combination of both non-universalities, 
scalars and gauginos, gives also rise to other interesting possibilities. In particular, 
these are choices giving rise to very light neutralinos. These arc obtained for exam- 
ple using the scalar non-universalities given in Eq. (2.2), and adding gaugino non- 
universalities. A representative case is 3 = 3 

As mentioned above, one can find in [19, 21], and [27], a detailed study of all these 
cases, and in particular of the computation of ncutralino-proton cross sections o'^o_„ 
and gamma-ray fluxes $-y, respectively. As an example, we show in Fig. 1 the values 
of c^o.p as a function of the neutralino mass ?w^o, for the case a) with tan/3 = 35, 
A = and /Lt > 0. The areas that will be covered in the future by two very sensitive 
detectors, EDELWEISS II and GENIUS, arc shown. Likewise, for the same case, we 
show in Fig. 2 the values of as a function of m^o^. In addition to EGRET, we show 
in the figure the region that will be covered by GLAST experiment. For the sake of 
definiteness, we use in our analysis of the gamma-ray flux an angular acceptance of the 
detector Ail = 10~^ sr, and a NFW profile [36]. The resuls can easily be extended to 
other profiles. For example, the flux produced by the Moore et al. profile [37] is about 
two orders of magnitude larger than the one produced by the NFW, and this is about 
two orders of magnitude larger than the flux produced by the Isothermal and Kravtsov 
et al. [38] profiles (see e.g. Table 1 in [27]). 

constructions in the type I string [32]. An analysis of direct detection of neutralinos in these scenarios was 
carried out in [33] . An analysis of the gamma-ray fluxes might also be interesting [34] . On the other hand, 
AMSB scenarios in an effective heterotic framework were studied in Ref. [25, 35]. 
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Let us remark that in the figures, only points allowed by all experimental constraints 
are shown. In particular, the lower bounds on the masses of the SUSY particles and on 
the lightest Higgs have been implemented, as well as the experimental bounds on the 
branching ratio of the b —>■ process and on a^^^^. Concerning the latter, we have 
taken into account the recent experimental result for the muon anomalous magnetic 
moment [39], as well as the most recent theoretical evaluations of the Standard Model 
contributions [40]. It is found that when e"'"e~ data are used the experimental excess 
in (51^ — 2) would constrain a possible SUSY contribution to be a^^^^ = (27.1 it 10) x 
10-^°. At 2(7 level this implies 7.1 x 10"^° < a^^SY < i ^ IQ-io, and therefore only 
cases with ^ > are interesting. It is worth noticing here that when tau data arc used 
a smaller discrepancy with the experimental measurement is found. In order not to 
exclude the latter possibility we will show the line a^^^^ = 7.1 x lO"^'^ in the figures 
throughout the paper. The region of the parameter space to the right of this line will 
be forbidden (allowed) if one consider electron (tau) data. 

Concerning the astrophysical bounds on the relic density, we show in the figures 
with black stars points with 0.094 < ^^oh'^ < 0.129, and therefore favoured by WMAP. 
Points with 0.129 < Q^oh"^ < 0.3 are shown with magenta triangles. We also consider 
the possibility that not all the dark matter is made of neutralinos, allowing 0.03 < 
n^oh"^ < 0.094. In this case we have to rescale the density of neutralinos in the galaxy 
p{r) by a factor Q^oh^ / 0.094. Points corresponding to this possibility are shown with 
blue boxes. For details on how the rest of the experimental bounds are implemented, 
and which programs are used to carry out the evaluation, see [27]. 
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Figure 1: Scatter plot of the scalar neutralino- proton cross section function of the neutralino 

mass m-^o in the non-universal scalar case a) Si — 0, S2 — 1, discussed in Eq. (2.2), for tan/3 — 35 and 
A = 0, with /i > 0. The region to the right of the dashed line corresponds to a^^^^ < 7.1 x 10^^", and 
would be excluded (allowed) by electron (tau) data. The points in the figure fulfil all other experimental 
constraints and they have 0.129 < ilj^oh^ < 0.3 (magenta triangles), 0.094 < ft^oh'^ < 0.129 (black stars), 
and 0.03 < ft^oh'^ < 0.094 (blue boxes) with the appropriate rescaling of the density of neutralinos in the 
galaxy as discussed in the text. ^From top to bottom, the upper area bounded by the first (second) solid 
line will be anahzed by the projected EDELWEISS 11 (GENIUS) experiment. 




Figure 2; The same as in Fig. 1 but for gamma-ray integrated flux for a threshold of 1 GeV versus the 
neutralino mass rn^^o . A NEW profile is used. Now, from top to bottom, the first solid line corresponds to 
the signal reported by EGRET, and the upper area bounded by the second solid line will be analyzed by 
GLAST experiment. 
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3 A comparison between direct and indirect de- 
tection 



As discussed in detail in previous sections there are many experiments in preparation 
for the direct and indirect detection of dark matter, and they wih be able to cover 
an important range of the parameter space of SUSY models. Thus, in principle, it 
might be possible that this means of detection will be able to probe a dark matter 
signal before the end of the decade. Therefore, it becomes interesting to compare 
both kinds of experiments. In particular, we will compare the direct detection of 
neutralinos, through their elastic scattering on target nuclei through nuclear recoils, 
with the indirect detection of neutralinos, through the gamma rays produced by their 
annihilation in the galactic center. 

Let us start with the simple mSUGRA scenario where the soft terms are universal 
at the GUT scale. In Fig. 3 we show in red the points fulfilling all experimental and 
astrophysical constraints discussed in the previous section, in the parameter space (m, 
M). The black colour on top of these indicates those points that will be testable by 
the experiments. In particular, on the left and central frames we use two very sensitive 
experiments for the direct detection of dark matter, EDELWEISS II and GENIUS, 
respectively. On the right frame we use the indirect detection experiment GLAST, 
assuming the NFW density profile for the theoretical computation of the signal. This 
figure is for tan /3 = 35, and throughout the paper we are using ^ = and > 0. 

Although indirect detection experiments will be able to test larger regions of the 
parameter space than direct detection ones, these regions are in any case small, and 
most of the points will avoid detection. Let us recall, however, that we are working 
with a NFW profile in our computation of gamma-ray fluxes. Since the Moore et al. 
profile is about two orders of magnitude larger, one should expect more detectability. 
This is summarized for different density profiles in Fig. 4, where Log( g^^g^^^y ) as a 
function of m-o is shown. For direct detection the signal consists of (t-o_„ and the 
sensitivity used in the plot is the one of the GENIUS experiment (see Fig. 1). For 
indirect detection the signal consists of $-y and the sensitivity used is the one of the 
GLAST experiment (see Fig. 2). The solid line separates points which are detectable 
(upper region) from those which are not (lower region). 

We have also analyzed very large values of tan /3, such as 50. In this case the differ- 
ence of sensitivity between direct and indirect detection experiments is more obvious. 
Whereas the region of the parameter space that will be tested by GENIUS is very 
similar to the one of the previous case with tan/3 = 35, the region tested by GLAST 
increases significantly Note, finally, that EDELWEISS II will not be able to detect any 
neutralino in the framework of mSUGRA. 

On the other hand, for non-universal Higgs masses the experimental situation 
changes significantly. Larger regions of the parameter space become available, mainly 
for GENIUS and GLAST. This can be seen in Figs. 5 to 10, which can be compared 
with Figs. 3 and 4 of mSUGRA. This effect is clearly stronger for the indirect detection. 
For example, only small regions of the parameter space will be tested by GENIUS for 
cases b) and c), whereas most of parameter space will be tested by GLAST. 

The situation for non-universal gauginos is also qualitatively similar, and we show 
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this in Figs. 11 and 12 

Finally, the other three representative cases mentioned in the previous section, with 
a very light neutralino mass, can be seen in Figs. 13 to 18. The same comments as 
above, for non-universal scalars, can be applied here. 
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Figure 3: Points fulfilling all experimental constraints with 0.1 < fl^oh'^ < 0.3 (red) in the parameter 
space of the mSUGRA scenario (m, M) for tan/3 = 35, A ^ and /i > 0. The black colour on top of these 
indicates those points that will be testable by EDELWEISS II (left), GENIUS (center), and GLAST (right) 
assuming the NFW density profile in the theoretical compution. The region to the right of the dashed line 
corresponds to a^^^^ < 7.1 x 10^ ^''j and would be excluded (allowed) by electron (tau) data. 
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Figure 4: Log( g^nsitivUy ) ^ function of the neutralino mass m^o for points fulfilling all experimental 
constraints with 0.094 < 17-^o/i^ < 0.129. The solid line separates points which are detectable (upper region) 
from those which are not (lower region). Black stars are associated to the case of direct detection, where 
the signal corresponds to cr^i-p ^^'^ sensitivity used is the one of GENIUS. The rest of the points are 
associated to the case of indirect detection through gamma-rays, where the signal corresponds to $.y and 
the sensitivity used is the one of GLAST. Results for different density profiles are shown: Moore et al. (blue 
boxes), NFW (red diamond), Kavtsov et al. (green triangle). For each profile, the region to the right of the 
dashed hue corresponds to a^^^^ < 7.1 x 10^^", and would be excluded (allowed) by electron (tau) data. 
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Figure 5: The same as in Fig. 3 but for the non- universal scalar case a) 6i ~ 0, 62 ~ I, discussed in 
Eq. (2.2). 
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Figure 6: The same as in Fig. 4 but for the non- universal scalar case a) Si — 0, 82 — I, discussed in 
Eq. (2.2). 
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Figure 7: The same as in Fig. 3 but for the non-universal scalar case b) Si — —1, 62 = 0, discussed in 
Eq. (2.2). 
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Fi gure 8: The same as in Fig. 4 but for the non-imiversal scalar case b) Si — —1, S2 = 0, discussed in 
Eq. (2.2). 
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Figure 9: The same as in Fig. 3 but for the non- universal scalar case c) 8\ = —1, ^2 = 1, discussed in 
Eq. (2.2). 
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Figure 10: The same as in Fig. 4 but for the non-universal scalar case c) 5i 
Eq. (2.2). 
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Figure 11: The same as in Fig. 3 but for the non-universal gaugino case d) 62 = 0, ^3 = —0.4, discussed in 
Eq. (2.4). 
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Figure 12: The same as in Fig. 4 but for the non-universal gaugino case d) S2 = 0, S'^ ~ —0.4, discussed in 
Eq. (2.4). 
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Figure 13: The same as in Fig. 3 but for the non-universal scalar and gaugino case Si ~ 0, S2 ~ 1, 
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Figure 14: The same as in Fig. 4 but for the non-universal scalar and gaugino case Si — 0, 52 — 1, 
5'. = S'^ 3. 
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Figure 15: The same as in Fig. 3 but for the non-universal scalar and gaugino case Si — —1, 62 — 0, 
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Figure 16: The same as in Fig. 4 but for the non-universal scalar and gaugino case Si = —1, 82 = 0, 
S'=S'.^ 3. 
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Figure 17: The same as in Fig. 3 but for the non- universal scalar and gaugino case Si — —1, 62 — 1, 
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Figure 18: The same as in Fig. 4 but for the non-universal scalar and gaugino case Si 
S'=S'.^ 3. 
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4 Conclusions 



There are two kinds of experiments in order to detect dark matter. Direct detection 
experiments through the elastic scattering of dark matter particles on target nuclei 
through nuclear recoils, and indirect detection experiments through the analysis of 
dark matter annihilation products. We have compared both kinds of experiments for 
neutralino dark matter, paying special attention in the indirect detection case to the 
gamma-ray fluxes produced by the annihilation of neutralinos in the galactic center. 

The discussion has been carried out using the parameter space of a general SUGRA 
theory, where not only universal soft scalar and gaugino masses are allowed, but also 
non-universalities may be present in both sectors, increasing significantly the theoreti- 
cal cross sections. For the comparison we pay special attention to points which will be 
testable by very sensitive experiments, such as e.g. EDELWEISS II and GENIUS for 
direct detection, and GLAST for indirect detection. 

In particular, we have determined the regions of the parameter space that will be 
testable by those two kinds of experiments and compare between themselves. We can 
summarize the results as follows. Interesting direct detection experiments, such as e.g. 
EDELWEISS II, will only be able to cover a extremely small region of the parameter 
space. More powerful detectors such as GENIUS will cover larger regions, but still 
not too large compared with the parameter space of general SUGRA. On the contrary, 
indirect detection experiments will be able to test important regions of the parameter 
space of SUGRA. 

Of course, a signal in a direct detection experiment would be a definitive proof of 
the existence of dark matter. However, we should keep in mind that this might be 
difficult, depending on the region of the parameter space where the neutralino dark 
matter is attached. Therefore, at the end of the day, indirect detection experiments 
might play a crucial role. 
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